Room-temperature exciton-polariton superfluid quantum interference device and quatron-polariton superconducting quantum interference device

ABSTRACT

The present disclosure relates to Superfluid QUantum Interference Devices (SQUIDs) that measure phase differences existing in quasi-particles or matter-wave systems, and the related techniques for their use at room-temperatures. These Bose-Einstein Condensation interferometry techniques include quantum scale metrology devices such as quasi-particle based linear accelerometers, gyroscopes, and Inertial Measurement Units that incorporate such interferometers. In the presence of additive white Gaussian noise, estimates are made for the Bias Instability, Angle Random Walk, and Velocity Random Walk of the device for purposes of quantum inertial sensing. Moreover, this disclosure relates to SQUIDs based on charged quasi-particles that can, in turn, be used to construct quantum computing elements such as quantum transistors, and quasi-particle circuits at room-temperatures. These quasi-particle circuits can be used to build analogs of electronic circuit elements, and offer an alternative to traditional electronics. Using a quasi-particle circuit, hysteresis can be achieved and controlled to build these new devices.

I. CROSS REFERENCE TO RELATED APPLICATIONS

This application claims priority to U.S. Provisional Application No. 62/275,953, filed Jan. 7, 2016, U.S. Provisional Application No. 62/343,629, filed May 31, 2016, U.S. patent application Ser. No. 15/402,149, filed Jan. 9, 2017, and U.S. patent application Ser. No. 16/358,071, filed Mar. 19, 2019, the contents of these applications are expressly incorporated herein in their entirety.

II. FIELD

The present disclosure relates to quantum interference devices that measure phase differences (for example) existing in quasi-particles and the related techniques for their use at room-temperatures.

III. BACKGROUND

In order to take measurements or extract phase information, interferometry relies on the superposition (i.e., interference of waves) that are typically electromagnetic. Interferometry has important applications in a wide range of scientific and engineering pursuits from astronomy to quantum mechanics, and particle physics. Generally speaking, interferometry makes use of the superposition of waves in such a way as their combination has some meaningful property that can be measured by way of the interference pattern created, which is indicative of the phase difference between the two or more waves. Waves that are in-phase undergo constructive interference while waves that are out of phase will undergo destructive interference. Waves that are neither completely in-phase, nor completely out-of-phase will create an interference pattern, which can be used to determine the relative phase difference between the waves. These phase differences can become apparent in the macroscopic phase of quantum condensed matter systems, thereby making them of interest for applications as quantum interference devices.

IV. BRIEF DESCRIPTIONS OF THE DRAWINGS

Implementations of the present technology will now be described, by way of example only, with reference to the enclosed figures, wherein:

FIG. 1 is an isometric view of an example of a polariton Superfluid Quantum Interference Device (pSQUID) and a diagrammatic view of an experimental configuration for producing vortex-antivortex superpositions in accordance with the present disclosure;

FIG. 2 is an isometric view of an example of a quatron Superconducting QUantum Interference Device (qSQUID) in accordance with the present disclosure;

FIG. 3 is a time evolution of a polariton superfluid current in a ring geometry in the presence of a Josephson junction in accordance with the present disclosure;

FIG. 4 is a diagrammatic view of a time evolution of a polariton Bose-Einstein condensate (BEC) in a disordered potential landscape in accordance with the present disclosure;

FIG. 5 is a diagrammatic view of a time evolution of a polariton BEC in a ring potential landscape in accordance with the present disclosure;

FIG. 6 is a diagrammatic view of a time evolution of a quatron BEC in a ring potential landscape with Josephson junctions in accordance with the present disclosure;

FIG. 7 is an Allan variance curve corresponding to the polariton accelerometer in the presence of additive white gaussian noise (AWGN) in accordance with the present disclosure;

FIG. 8 is an Allan variance curve corresponding to the polariton gyroscope in the presence of additive white Gaussian noise (AWGN) in accordance with the present disclosure;

FIG. 9 is a conceptual description of a polariton Inertial Measurement Unit (pIMU) in accordance with the present disclosure.

V. DETAILED DESCRIPTION

The type of devices discussed herein will permit practical implementations of quantum memory, digital noise filters (e.g., Schmitt triggers), magnetometers (e.g., Superconducting QUantum Interference Devices or SQUIDs), qubits, qudits, and quantum Inertial Measurement Units (IMUs). For example, ultra-cold atomic gases offer highly attractive physical systems for realizing such quantum devices, as they offer an exquisite experimental control that originates from quantum optics, such as vortex or soliton creation in a BEC matter wave. However, one of the drawbacks of ultracold atom technology is that in order to achieve BEC, conditions need to be cooled to nano-kelvin temperatures, which can prove quite difficult and expensive. In order to overcome this challenge, the present disclosure employs quasi-particles, which can form BEC matter waves at room temperatures. Examples of the quasi-particles are exciton-polaritons, and quatron-polaritons. Exciton-polaritons are half-light, half-matter quasi-particles that result from combing quantum well excitons (i.e., electron-hole pairs) and photons via planar. Fabry-Perot resonators, distributed Bragg reflector (DBR) microcavities, or photonic crystals, for example. Negatively charged quatrons are composed of a hole and three electrons in a semiconductor. Because quatrum-polaritons (charge=−2e) and exciton-polaritons (charge=0) are optically accessible bosons, they offer an entirely new arena for engineering BECs at room temperature via quantum optics and a new configuration for achieving novel quantum devices such as optoelectronics (e.g., quantum circuits, ultra-fast optical switches), power production (e.g., hybrid organic-inorganic solar cells), and medicine (e.g., compact terahertz lasers). Herein, the disclosure utilizes the superfluid property of polariton condensates, and the superconducting property of quatron condensates, for achieving these examples according to the present disclosure.

It will be appreciated that for simplicity and clarity of illustration, where appropriate, reference numerals may have been repeated among the different figures to indicate corresponding or analogous elements. In addition, numerous specific details are set forth in order to provide a thorough understanding of the examples described herein. However, it will be understood by those of ordinary skill in the art that the examples described herein can be practiced without these specific details. In other instances, methods, procedures and components have not been described in detail so as not to obscure the related relevant feature being described. The drawings are not necessarily to scale and the proportions of certain parts may be exaggerated to better illustrate details and features. The description is not to be considered as limiting the scope of the examples described herein.

Phase difference that is measured via interferometry is typically created by a difference in path length between a sample and reference beam that are created by splitting a single source of light, such as a laser. Even small differences in the path length of the sample beam can result in a phase difference in the detected, recombined beam.

Phase differences created by path length differentials are not the only phase change that can be measured by an interferometer. For example, a phase change can also result from a difference in refractive index between two paths of the same length.

Additionally, interferometers can also measure the phase differences created by differences in angular momentum existing in quantum vortex-antivortex superposition states, for example. The Orbital Angular Momentum (OAM) of light results in a so-called helical mode, which can be used to create these quantum vortex-antivortex superposition states in BECs. OAM states can occur naturally in light; however, there are a number of methods for creating OAM states in the laboratory. One method for creating an OAM state employs a Spiral Phase Plate (SPP). SPPs can be made of plastic or glass, for example, and are constructed such that the thickness of the material increases in a spiral pattern in order to imprint a phase gradient on light passing through it, which can in turn be transferred to a BEC. Alternatively, these OAM vortex-antivortex superposition states can be obtained via metastable condensation, where the condensate is placed in a circular narrow channel, i.e., ring potential. For the metastable condensation technique, the electrical injection of the condensate removes the need for optics. By eliminating the vortex core with the use of a ring potential, high OAM states become ultrastable, which is beneficial from a sensitivity point of view. Moreover, Josephson junctions can be placed in these ring potentials, thereby constituting quasi-particle analogues of SQUIDs. The magnitude of the barrier height of the junction creates a quasi-particle analogue of a Superconductor-NormalMetal-Superconductor (SNS) junction.

The Gouy phase and Sagnac phase are often measured in interferometer systems. The Gouy phase is the phase that a wave acquires as it changes its radius of curvature via diffraction. The Gouy phase of optical beams can be observed by the interference patterns that are created, i.e., superimposing a fundamental Gaussian beam going through a focal point with a plane wave. The Gouy phase and Sagnac phase can also originate from the interference of two counter-rotating OAM states.

Sagnac interferometers are often used for measuring angular velocity. The Sagnac interferometer splits the light from a single collimated source by way of a partially reflective mirror (e.g., beam splitter) into two beams that follow the same path but in the opposite directions and can be thought therefore to enclose an area. The relative phases of the two exiting beams, and thus the position of the interference pattern are shifted according to the angular velocity of the beams. This is also referred to as the Sagnac Effect.

Such interferometer arrangements have been incorporated into laser gyroscopes, and fiber optic gyroscopes, and the relative sensitivity of such gyroscopes arises from the invariance of the speed of light in all inertial frames of reference.

The Mach-Zehnder interferometer (MZI) is a type of Sagnac interferometer. An MZI also splits the light from a single collimated source using a partially reflective mirror (e.g., beam splitter). The two resulting beams are sometimes referred to as the “sample beam” and the “reference beam,” as discussed above. Often for experimental purposes a sample (that alters the phase of this beam) is placed within the path of the sample beam.

Each sample and reference beam is reflected by a mirror; the two beams then pass through a second partially reflective mirror (e.g., Beam Splitter or BS) to enter one or two detectors. Depending on the presence (and nature) of a sample within the sample beam, different interference patterns are created at the detectors.

MZIs can be used for precision phase measurement. The sensitivity of MZIs for phase measurement is limited by the shot noise of the photons which scales as the inverse of the square root of the number of photons, although the shot noise can be suppressed by using squeezed states of light at the input of the MZI. Alternatively, the signal can be enhanced by optical parametric amplifiers. Combined with the squeezed state input, the sensitivity can be improved further due to the noise reduction, thereby permitting the “Heisenberg limit,” which scales as the inverse of the number of photons.

Optical interferometers depend, to some extent, on the wavelength of the optical beam used, and thus are limited in sensitivity. Moreover, for example, optical Sagnac interferometers depend on the area of the configuration meaning that the sensitivity to rotation is proportional to the area circumscribed by the counter-rotating optical beams. Consequently, manufacturing an optical Sagnac interferometer in a compact form, thus having a lower relative circumscribed area, can come at the cost of a decrease in measurement sensitivity for a given wavelength. However, one advantage of optical interferometers is that they operate at room temperatures. Furthermore, such optical interferometers are relatively inexpensive.

An alternative to optical interferometers are atomic-based interferometers. Atomic interferometers make use of the properties of ultra-cold atomic gases and offer highly attractive physical systems for realizing very sensitive (i.e., m_(eff)c²/

ω˜10¹⁰), quantum-scale metrological devices for applications such as gyroscopes that employ the Sagnac effect. Atomic systems offer exquisite control that originates from quantum optics, such as vortex creation in a BEC matter wave, and are highly coherent, meaning the interference visibility can be used for interferometry.

Atomic BEC was first achieved in 1995 and inspired a massive worldwide research effort. Today, atomic BECs (for example, based on Helium isotopes) are being considered for applications in a variety of fields, such as quantum simulation, quantum information, as well as quantum metrology. BECs can exhibit a number of unusual physical properties, such as superfluidity, vortices, and attractive interactions.

The rotational measurement sensitivity of atomic matter wave interferometers per unit area typically far exceeds that of optical interferometers. Unfortunately, achieving BEC in atomic gases is very difficult and can typically only occur at nano-kelvin temperatures, rather than at room temperatures (i.e., approximately 300 Kelvin). As a result, atomic BEC interferometers can be impractical and extremely expensive.

Atomic BECs are not the only type of BECs that can be created. BECs can also be formed of quasi-particles. Quasi-particles are so-called emergent phenomena that occur when a system in a solid behaves as if it contained different weakly interacting particles in free space. For example, as an electron travels through a semiconductor, its motion is disturbed by its interactions with all of the other electrons and nuclei; however, it approximately behaves like an electron with a different mass (i.e., effective mass) travelling unperturbed through free space. Such an electron can be referred to as an electron quasi-particle.

As another example, the aggregate motion of electrons in the valence band of a semiconductor can be modeled as if the semiconductor instead contained positively charged quasi-particles called electron-holes or holes.

Polaritons are a class of quasi-particles resulting from the strong coupling of electromagnetic waves (e.g., light) within an electric or magnetic dipole carrying excitation. There are several types of polaritons. Exciton-polaritons result from the coupling of visible light with an exciton, i.e., an electron-hole pair.

Now with specific reference to FIG. 1, FIG. 1 illustrates a system 100 utilizing a MZI configured for producing vortex-antivortex superpositions in polariton BECs according to one aspect of the current disclosure. A laser 101 generates a Transverse Electromagnetic (

_(=0,p=0)) Gaussian beam profile. This light is passed through a SPP 102, generating an OAM state of light. Once the OAM state of light has traveled through the MZI with a Dove prism (DP) 103 and mirrors (M) 104 and M 105, an OAM superposition state of light is made. The resulting OAM superposition state of light can be used to generate an interference pattern in a polariton BEC matter wave within the DBR microcavity 106. This interference pattern is imaged by a detector such as a charge-coupled device (CCD) 107, and can then be used to determine the Sagnac or Gouy phase of the polariton BEC, thereby constituting a pSQUID. Here, it should be pointed out that an annular potential can be placed in the DBR microcavity 106 such that a polariton superfluid ring can be generated. As such, Josephson junctions can also be placed in the annular potential in the DBR microcavity 106.

As seen in FIG. 2, negatively charged quatron-polaritons are composed of a hole and three electrons in a triple-layered semiconductor microcavity, where a two-dimensional electron gas (2DEG) 205 layer is sandwiched between two BEC layers 205, for example. If coupled to superposition states of light from an OAM emitter 201, for example, they become stable. Because quatron-polaritons are optically accessible bosons, they offer a way for the OAM of light to be coupled to a BEC (e.g., in semiconductors), and are robust at room temperature. Moreover, this facilitates the entanglement of quasi-particle BECs, as OAM states of light can be entangled.

Hereinafter, “exciton-polaritons” will be referred to as “polaritons” for brevity and simplicity, while “quatron-polaritons” will be referred to as “quatrons” for similar reasons.

Polariton BECs are formed in semiconductor microcavities, planar dielectric Fabry-Pérot optical cavities, or photonic crystals, for example, and can be described as half-matter, half-light excitations that are typically created using well-known laser injection methods.

Due to the de Broglie wavelength of the polaritons, a polariton superfluid is a type of BEC that can be achieved at much higher temperatures than atomic BECs. Superfluidity is the resistanceless flow of a fluid, and is fundamental to describing many unusual physical effects. Such effects include the superfluid fountain, the ability of a superfluid to resist rotation while inside a rotating container, and the ability to “defy” gravity.

Polariton BECs are a non-equilibrium phenomenon that can be achieved at room temperatures for suitable materials that support polaritons. In materials where polaritons can be operated at room temperature, interactions tend to be weaker; however more recently, the nonlinearities have been observed to be significant. Despite the non-equilibrium nature of the polariton BEC, it shares many of the properties of an atomic BEC, such as superfluidity, stable vortex formation, and high spatial coherence. Polaritons can have very high critical temperatures clue to their exceedingly low effective mass (i.e., 10⁻⁵ times the bare electron mass). As discussed below, polariton BECs can be used to construct pSQUIDs.

Quatron BECs, as described above, carry negative charges and are useful because they exhibit superconducting properties. Again, quatrons are quasi-particles composed of a hole and three electrons which form a stable BEC when coupled to light in triple quantum layer structures in semiconductor microcavities, as seen in FIG. 2. This allows for both the greater experimental control found in quantum optics, and the benefits of matter wave systems, such as superconductivity and coherence. Moreover, due to the extremely small effective mass of the quasi-particles, quatrons can be used to achieve superconducting BECs at room temperature. As discussed below, quatron BECs can be used to construct qSQUIDs.

Hysteresis is a general feature of systems where the energy has two (or more) local minima separated by an energy barrier (e.g., a Josephson junction). For ring-shaped, quasi-particle BECs, these minima represent stable flow states of the system (e.g., counter-rotating electron-hole currents), and their energies can be altered via rotating the BEC with an angular velocity Ω, or rotating the Josephson junctions. As the quasi-particles interact with one another, an energy barrier can appear, creating multiple flow states. The energy barrier stabilizes the flow, making the quasi-particle BEC a superfluid, or superconductor.

These ring-shaped superfluids (or superconductors) exhibit both hysteresis and a critical angular velocity. As the rotation of the ring generates a quantized flux, the present device has applications as a novel flux qubit, for example. Additionally, owing to the optical coupling, these ring-shaped quasi-particle superfluids (or superconductors) can be entangled.

Superfluidity in polariton condensates are analogous to other systems (e.g., superfluid helium or atomic condensates) which sets them apart from other solid-state quantum systems. It has recently been demonstrated that hysteresis exists between quantized circulation states in a ring of superfluid atomic BEC obstructed by a rotating weak link. This was accomplished with a tunneling barrier such that there exists a region of low density in the superfluid ring.

Thus, according to one or more examples of the present disclosure, polariton-based interferometry can be conducted through the use of room-temperature polariton BECs. More particularly, OAM superpositions of light can be created and used as a seed to imprint these OAM state superpositions onto room-temperature polariton BEC clouds. The Sagnac or Gouy phase, for example, of the superposition of OAM states can be detected via the resulting interference pattern.

Through novel interferometer arrangements relying on the polariton BEC superposition states, sensitivity can be achieved in compact form factors exceeding what is possible with optical or atomic interferometers of the same occupied area. Moreover, because of the low relative effective mass of polaritons, the quasi-particle BEC can be created and sustained at room temperature conditions, unlike atomic BEC interferometers. As a result, such devices can achieve exquisite sensitivity while being chip-scale manufacturable and practical for device miniaturization purposes.

These polariton interferometers, with modifications discussed as part of other examples of the present disclosure, can be constructed to operate as linear accelerometers and gyroscopes, and thence combined to construct IMUs and other devices for related applications.

Other examples of the present disclosure permit the construction of quatron-based quantum transistors for use as elements in quantum-based logic and computations. The charged supercurrents found in room-temperature quatron BECs can be used to replace the voltage state representations of classical transistors, as discussed below, for example. In classical electronic transistors, information is encoded and accessed by manipulating and reading voltages on transistor switches located inside integrated circuits. Each electronic transistor is set to a voltage state of either 0 or 1 by a bus to encode bits of information. Devices based on qSQUIDs, however, can be used to construct quantum transistors, logic blocks, and quantum computing elements as they are set to a voltage states 0 and 1 simultaneously. This allows qSQUIDs to function as quantum bits (i.e., qubits), or qudits, as they encode the quantum information from the OAM states of light in the quatron BEC.

According to some examples of the present disclosure, the type of devices described will permit the construction of quantum memory components, digital noise filters (such as Schmitt triggers), magnetometers (e.g., SQUIDs), qubits, qudits, quantum IMUs, quantum Inertial Navigation Systems (INSs), and different varieties of interferometers that are based on quasi-particle physics rather than the well-known optical or atomic counterparts. Moreover, the examples of the present disclosure can be operational at room temperature (i.e., approximately 300 Kelvin). More particularly, OAM superpositions of light or electrical injection can be used as a seed to imprint superpositions states onto room temperature polariton BEC clouds. The Sagnac phase, Gouy phase, or magnetic flux quantum, for example, of the superposition states can then be detected via the resulting quantum interference patterns, or magnetic flux quantum.

According to other examples of the present disclosure, ring geometries are used where the polariton condensate is placed in a circular narrow channel potential well. An alternative method of forming the superposition states in a polariton BEC is via metastable condensation, where the ability of electrical injection of the BEC removes the need for optical injection.

These interferometers with modifications discussed as part of some examples of the present disclosure can be constructed to operate as linear accelerometers, gyroscopes, and hence combined to construct IMUs. Other examples of the present disclosure permit the construction of quatron based transistors for use as an element in quantum-based logic and computation.

Many lasers emit light coherent states which approximate a Gaussian profile with

=p=0. Using the TEM₀₀ fundamental transverse mode of the optical resonator and channeling it through a SPP, for example, OAM states of light can be generated. The resulting OAM state of light is sometimes referred to as a Laguerre-Gauss (

_(p)) mode. The

_(p) optical state with a topological charge

is

$\begin{matrix} {{\left. ❘ \right\rangle = {\left( \frac{2}{\pi{{\ell }!}} \right)^{1\text{/}2}\frac{1}{w}\left( \frac{r\sqrt{2}}{w} \right)^{\ell }{\exp\left( \frac{- r^{2}}{w^{2}} \right)}{\exp\left( {{- i}\;\ell\;\phi} \right)}{\exp\left( {i\left\lbrack {{kz} - {{kr}^{2}\text{/}\left( {2R} \right)}} \right\rbrack} \right)}{\exp\left( {i\;\phi} \right)}}},} & (1) \end{matrix}$ where r, ϕ, z are cylindrical coordinates, w is the beam waist, k is the wave number, R is the radius of curvature of the wave-front. Similarly, the |−

state can be obtained by replacing

with −

in Eq. (1). Although these states can be made with holographs or cylindrical lens mode converters and other means, SPPs offer a simple and convenient way to obtain a pure OAM state.

Once pure OAM states are created, they can be made into OAM superposition states. One way to accomplish this superposition is by using a MZI, where a DP 103 is used to transform the OAM state. By altering the reflectivity of the first BS 108 and second BS 109 of the MZI, or altering the phase 111, an arbitrary superposition of the OAM states can be obtained.

The generated superposition of OAM states can be used to coherently excite the polaritons in, for example, a DBR structure 106. DBRs are semiconductor devices used to facilitate polariton BECs. As an alternative to semiconductors, photonic crystals, or polymers are materials that will also support the polariton BEC, for example. By exciting the polaritons accordingly, the seed of the polariton condensate is reinforced by the standard process of exciton-polariton cooling and scattering. Once the characteristic interference of the OAM superposition has been transferred to the polariton BEC, the Sagnac phase can then be determined from the angular velocity of the BEC cloud, and the Gouy phase can be determined from the linear velocity of the BEC cloud at the CCD 107.

The source of coherent light is a TEM₀₀ emitting laser 101. The beam of laser light is channeled through the SPP 102 to create an OAM state of light or a LG

_(,p). This OAM state of light then travels through a first beam splitter (BS) 108. The partially reflective surface of BS 108 divides the beam into two arms according to its reflectivity.

The lower arm of the divided source of OAM state light passes through the DP 103. The DP 103 has the effect of creating a counter-rotating OAM state of light that serves as the antivortex to the OAM state of light exhibited by the upper arm of the MZI. The lower arm reflects from mirror (M) 104 and then is channeled into the second beam splitter BS 109, which in one example can be 50% reflective, 50% transmissive.

The upper arm reflects from mirror (M) 105 and is also channeled into the second beam splitter BS 109. In certain configurations, a sample can be placed at item 111, which may introduce a phase shift in the sample beam. The imaginary output 110 emerging from the second BS 109 is typically disregarded. The real output of the second BS 109 represents the superposition of the two OAM states. This superposition is channeled into the DBR microcavity 106, where (as discussed above) it serves as a seed polariton condensate.

The DBR 106 is itself pumped by a coherent light source 112. The output of source (e.g., laser) 112 is reflected by mirror (M) 113 to provide the hot polaritons for condensation by the DBR 106 using known techniques. The DBR 106, as will be discussed in more detail below, emits a luminescence that includes the interference pattern of the counter-rotating vortex-antivortex superposition states. This output from the DBR 106, containing the interference pattern, can be detected at 107, which is for example a CCD, and thus the Sagnac and Gouy phase of the polariton BEC can be determined using standard techniques. Alternatively, as will be discussed below other devices or configurations may be employed rather than a detector to achieve different measurements and effects.

The laser injection methods discussed above can be used to resonantly excite the polariton BEC. The phase of the light is directly imprinted on the polariton BEC system, thus the injected polaritons simply follow the same phase relation as the injected light. In this regard, there is no difference in using the OAM superposition states of light to the polariton BEC as the polaritons inherit their vortex nature from the light. However, once the initial coherent excitation is induced, the OAM superposition can be turned off and the microcavity system is excited by conventional polariton pumping. This can be performed by an off-resonant pumping where the excitors have a much higher energy than the polaritons, or by pumping at a high in-plane momenta. Either method creates a large population of uncondensed polaritons from which the polariton BEC is replenished.

Notably, apart from the initial seed polariton BEC which is induced coherently, the remaining population is not pumped with the OAM superposition state. Hence in the second phase where the conventional polariton BEC pumping is used, any OAM superposition state is a genuine steady state of the polariton BEC and possesses its own dynamics. Vortex-antivortex superpositions have a more complex phase relation than standard vortex configurations, and have the possibility to spontaneously separate spatially, or relax to a zero momentum state.

The spatiotemporal evolution of the non-equilibrium polariton system is described by the open-dissipative Gross-Pitaevskii (dGP) equation, expressed as:

$\begin{matrix} {{{i\frac{\partial{\psi\left( {x,t} \right)}}{\partial t}} = {\left\lceil {{{- \frac{\hslash}{2m_{eff}}}{\nabla^{2}{+ {V_{ext}(x)}}}} + {{\mathcal{g}}{{\psi\left( {x,t} \right)}}^{2}} + {\frac{i}{2}\left( {{P(x)} - \gamma - {\eta{{\psi\left( {x,t} \right)}}^{2}}} \right)}} \right\rceil{\psi\left( {x,t} \right)}}},} & (2) \end{matrix}$ where ψ(x,t) is a complex valued condensate order parameter, g is the polariton-polariton interaction constant, m is the polariton mass, V_(ext)(x) is a spatially dependent trapping potential energy, P(x) is a spatially dependent pumping rate, γ is the polariton loss rate, and η is the gain saturation. Throughout this disclosure the time scale is in units of 2ma²/

, where a is the length scale set by the experiment (e.g., size of trapping potential or pump spot diameter). For the case of a Josephson junction in the ring geometry, the external potential V_(ext)(x) will have two components. The first component of this potential is the static, trapping potential V _(trap)(x,y)=V ₀ −V ₀ exp[−2(√{square root over (x ² +y ²)}−R)² /w ²],  (3) where V₀, R, and w are set to approximate our model of the circular narrow channel potential. The second component is the potential that generates the Josephson junction, and is taken to be a Gaussian along the azimuthal direction and a rectangular potential in the radial direction. At θ=0 the form of the Josephson junction potential V_(JJ) is the product of a Gaussian along y and a rectangle along x:

$\begin{matrix} {{{V_{JJ}\left( {x,y} \right)} = {{\exp\left\lbrack {{- 2}y^{2}\text{/}w_{L}^{2}} \right\rbrack}{V_{up}\left( {x,{{- L}\text{/}2},a} \right)}{V_{down}\left( {x,{L\text{/}2},a} \right)}}},{where}} & (4) \\ {{{V_{up}\left( {x,x_{up},a} \right)} = \frac{{\tanh\left\lbrack {a\left( {x - x_{up}} \right)} \right\rbrack} + 1}{2}},{and}} & (5) \\ {{V_{down}\left( {x,x_{down},a} \right)} = {\frac{{\tanh\left\lbrack {a\left( {x_{down} - x} \right)} \right\rbrack} + 1}{2}.}} & (6) \end{matrix}$

Here, Eqs. (5) and (6) approximate step functions (i.e., V_(up) and V_(down)). The quantities x_(up) and x_(down) denote the points where the steps up and down occur, respectively, and a is a parameter that controls the sharpness of the jump. The actual potential is V_(JJ)(x′, y′) where x′ and y′ are the rotated coordinates: x′(x,y,θ)=x cos(θ)+y sin(θ),  (7) and y′(x,y,θ)=−x sin(θ)+y cos(θ).  (8) FIG. 3 illustrates a basic implementation of the polariton superfluid currents in a ring geometry in the presence of a Josephson junction. Cases shown are (a) t=0; (b) t=5; (c) t=10. Parameters used in Eq. (2) are γ=η=g=1, and R=10 in Eq. (3). A spatially dependent pump P(x, y)=P₀|ψ(x,y)|², with P₀=2 is used. Timescales are in units of 2ma²/

where a is the lengthscale unit.

FIG. 4 illustrates a stable and coherent vortex-antivortex superposition. Moreover, the vortex-antivortex superposition has the characteristic phase dependence of the lobes being out of phase with each other by a factor of π. The system quickly reaches steady state and little change is seen beyond t=10. The stability of the counter-rotating vortex-antivortex superposition can be attributed to the same reason as for atomic BECs: both the vortex and the antivortex require a density defect due to the phase velocity diverging at the vortex core. Thus, it is energetically favorable for both to share the same vortex core, effectively pinning them in the same location. Time evolution of a polariton BEC initially seeded with a vortex-antivortex superposition for a disordered potential. FIG. 4 includes timeframes for (a) initial state t=0; (b) intermediate state t=5; and (c) steady-state t=10. Timescales are in units of

/meV. The following parameters were used in the simulations:

g=0.05,

γ=1.0, and

η=0.1, where energies are in meV.

Generally, higher sensitivities of quantum interference are achieved by the use of high momentum states due to their relatively shorter wavelengths, which allow for higher resolutions. A variant of the vortex-antivortex superposition state within the scope of the present disclosure is to use ring geometries, where the polariton BEC is placed in a circular narrow channel potential, as seen in FIG. 5. In a BEC, high angular momentum states are energetically disfavored and typically multi-quantized vortices break up into many

±1 vortices. This can be understood from the high momentum states that are near the vortex core which have a rapid phase variation. By eliminating the vortex core with the use of a ring potential, high angular momentum states become ultra-stable, which is beneficial from a sensitivity point of view. To create such a higher angular momenta superposition, a metastable condensation technique is utilized, which can be accomplished by electrical injection of the polariton BEC. FIG. 5 is an example of the interference pattern of higher angular momentum states in a ring potential. This interference pattern can then be used to determine the Sagnac or Gouy phase, effectively making the superfluid polariton BEC a quantum inertial sensor. In FIG. 5, counterpropagating momentum currents are shown in ring geometries for (a)

=±1, t=0; (b)

=±1, t=10; (c)

=±5, t=0; (b)

=±5, t=10. Parameters used in Eq. (2) are γ=η=V₀=g=1, and R=5 in Eq. (3). A spatially dependent pump P(x,y)=P₀|ψ

(r, ϕ)|², with P₀=2 is used. Timescales are in units of 2ma²/

where a is the lengthscale unit.

From the interference pattern of the polariton BEC system as seen in FIGS. 4 and 5, when the device rotates in time, the two counter-rotating matter waves pick up different phases. The macroscopic quantum mechanical phase ϕ(t) accumulated over a given time duration can be calculated using the prescription of the Sagnac effect, or Gouy effect.

Superfluid polariton condensates can be described by a Landau two fluid model, with a macroscopic wavefunction Ψ(x,t)=√{square root over (ρ_(s))}·exp(iϕ(t)),  (9) where ρ_(s) is the superfluid density and ϕ(t) is the macroscopic phase. The superfluid mass content I across a Josephson junction can be written as I=I ₀ sin(ϕ(t)),  (10) and is analagous to a superconducting dc-Josephson equation, which is also valid for BECs. The phase difference ϕ(t) propagates in time according to a Josephson-Anderson phase equation

$\begin{matrix} \begin{matrix} {{\overset{.}{\phi}(t)} = {- \frac{\Delta\mu}{\hslash}}} \\ {{= {- \frac{m_{eff}\left( {{\Delta\; P\text{/}\rho_{s}} - {s\;\Delta\; T}} \right)}{\hslash}}},} \end{matrix} & (11) \end{matrix}$ where m_(eff) is the effective mass of the constituent superfluid particles (˜10⁻⁵ the electron mass for polaritons), ρ_(s) is the superfluid density, s is the entropy per unit mass, ΔP and ΔT are the pressure and temperature differentials, respectively. For polariton BECs, a phase transition to superfluidity can occur below the transition temperature at a given density, or above the transition density at a given temperature. The critical transition temperature for a two-dimensional Bose gas is governed by the relation

$\begin{matrix} {{T_{c} = \frac{2{\pi\hslash}^{2}n_{2}}{1.645m_{eff}k_{B}}},} & (12) \end{matrix}$ where κ_(B) is the Boltzmann constant. For a polariton gas with a number density of n₃=3.5×10¹¹ [cm⁻³], the two-dimensional concentration is n₂=0.5×10⁸ [cm⁻²], which corresponds to room-temperature T_(c)=300 K. Near this transition temperature, oscillations are due either to Eq. (10) or to a linear current-phase relation with 2π phase slips. If a polariton superfluid loop contains two Josephson junctions in the annular potential as given by Eqs. (3)-(6), and a chemical potential Δμ is applied across them, the current across each junction will oscillate as I₁ sin(ϕ₁(t)) and I₂ sin(ϕ₂(t)), respectively. By integrating the superfluid loop,

$\begin{matrix} \begin{matrix} {{\oint{{{\nabla\phi} \cdot d}\; 1}} = {{\phi_{1}(t)} - {\phi_{2}(t)} + {\phi_{ext}(t)}}} \\ {{= {2\pi\; n}},} \end{matrix} & (13) \end{matrix}$ where ϕ_(ext(t)) is an external phase shift, i.e., the Sagnac phase shift for rotation, or the Gouy phase shift for acceleration. Hence the amplitude of the total current can be written

$\begin{matrix} {{{I(t)} = {\left( {I_{1} + I_{2}} \right)\sqrt{{\cos^{2}\left( \frac{\phi_{ext}(t)}{2} \right)} + {\xi^{2}\mspace{14mu}{\sin^{2}\left( \frac{\phi_{ext}(t)}{2} \right)}}}}},} & (14) \end{matrix}$ where ξ=(I₁−I₂)/(I₁+I₂) describes the asymmetry across the Josephson junction. The pSQUID is a charge neutral and room-temperature operational version of a de-SQUID. However, a charged supercurrent device can be obtained by implementing Josephson junctions with quatron-polaritons. When the pSQUID is placed in a rotating frame, the polariton superfluid rotates in a direction normal to the partition walls containing the Josephson junctions, thereby functioning as a quantum inertial sensor.

The Gouy phase is the phase that a wave acquires as it changes its radius of curvature due to diffraction. The Gouy phase of optical beams can be observed by interference, i.e., superposing a fundamental Gaussian beam going through a focal point with a plane wave. Hence,

$\begin{matrix} {{\phi_{Gouy}(t)} = {\left( {\alpha + 1} \right)\mspace{14mu}\tan^{- 1}\left\lceil \frac{z(t)}{z_{R}} \right\rceil}} & (15) \end{matrix}$ is the Gouy phase, z_(R) is the Raleigh range, and α=2p+|

|. When considering a superposition of two LG modes with topological charges

and −

and the total state is |ψ

=sin(θ)|

₁=

+cos(θ)|

₂=−

exp(iδ),  (16) where δ is the relative phase between the two beams. Eq. (16) describes the amplitude of the beam in FIGS. 1 and 2 after it has passed through the BS 109 or the OAM emitter 201. For a two-dimensional condensate such as a polariton BEC, the angular change in the vortex position between two points along the z-direction is

$\begin{matrix} {{{\overset{.}{\phi}}_{Gouy}(t)} = {\frac{{\ell_{2}} - {\ell_{1}}}{\ell_{2} - \ell_{1}}{\left\lceil \frac{\overset{.}{z}t}{z_{R}} \right\rceil.}}} & (17) \end{matrix}$

It has recently been shown that laser injection methods can be used to generate counter-rotating polariton BEC superpositions as seen in Eq. (16). Via Eq. (12), an exciton-polariton (i.e., polariton) superfluid is a type of BEC that can be achieved at much higher temperatures than typical BECs. Assuming p=0, the dynamic phase δ is constant in Eq. (16), and

₁=

,

₂=−

then

$\begin{matrix} {{{\overset{.}{\phi}}_{Gouy}(t)} = {\left( {\alpha + 1} \right){\left\lceil \frac{z_{R}\overset{.}{z}t}{z_{R}^{2} + z^{2}} \right\rceil.}}} & (18) \end{matrix}$ For a two-dimensional condensate, Eq. (18) is

$\begin{matrix} {{{\overset{.}{\phi}}_{Gouy}(t)} = {\left( {\alpha + 1} \right){\left\lceil \frac{\overset{.}{z}t}{z_{R}} \right\rceil.}}} & (19) \end{matrix}$ The signal-to-noise ratio (SNR) and the estimates for achievable sensitivity in a polariton Gouy accelerometer is discussed below. The SNR is given by the phase shift, i.e., Eq. (19), to the shot noise that depends on the number of photons N at the detector from the condensate per second,

$\begin{matrix} {{SNR} = {\frac{{\overset{.}{\phi}}_{Gouy}(t)}{\sqrt{1\text{/}N}}.}} & (20) \end{matrix}$ The final SNR has a maximum at a certain rate of photon detection, and the sensitivity ż_(min) is calculated by setting SNR=1. Thus, for the OAM superposition,

$\begin{matrix} {{\overset{.}{z}}_{\min} = {\frac{z_{R}}{t\sqrt{N}\left( {\alpha + 1} \right)}.}} & (21) \end{matrix}$

Typical polariton densities are in the region of 10⁹ [cm²], and a polariton spot size of 100 [μm²] gives a polariton number in the region of 10³. The lifetime of the polaritons is in the region of ˜10 [ps], hence the number of signal photons can be estimated as N˜10¹⁴ [s⁻¹] for polaritons. The Rayleigh range z_(R) is ˜1×10⁴ [cm]. This gives an estimate for the polariton vortex-antivortex superposition as ż _(min)˜5×10⁻⁹[m·sec⁻¹·Hz^(−1/2)].  (22)

Here, it should be pointed out that if |

₂|≠|

₁| in Eq. (17), Gouy phase changes can be detected via rotations of the quantum interference. However, for a vortex superposition where |

₂|=|

₁|, Eq. (19) can be used to determine the Gouy phase shift, as there will be no rotations of the quantum interference. OAM superposition states of light, i.e., Eq. (16) can be generated by passing a TEM₀₀ beam mode through a SPP, followed by a MZI with a DP. The SPP converts the TEM₀₀ to a LG

_(,p) beam mode Eq. (1) of quantized OAM with quantum numbers

and p, and the DP changes the winding of the quantized OAM from

to −

. In addition to generation using SPPs, cylindrical lens mode converters, spatial light modulators, and computer generated holograms are alternative methods to generate integer quantized OAM superposition states. Alternatively, an OAM emitter 201 can be used to generate the OAM superposition states of light, as seen in FIG. 2. The LG

_(,p) beam modes form a complete orthonormal basis set of solutions for paraxial light beams commonly found in lasers. These modes can be expressed by Eq. (1).

Typical polariton effective masses are m_(eff)˜10⁻⁴m_(e), where m_(e) is the electron mass. By comparing the rest mass energy of polaritons to the photon energy, m_(eff)c²/

ω˜1. In a polariton BEC where the phase around a vortex is fixed topologically by the OAM winding number

, a phase of 2π

accumulates about the vortex. Since the Sagnac phase originates from the interference of two counter-rotating orbital angular momentum, the Sagnac phase is written ϕ_(Sagnac)(t)=2

{dot over (ϕ)}t.  (23) For polariton BECs the Sagnac phase can be written

$\begin{matrix} {{{\phi_{Sagnac}(t)} = {{N(t)}\frac{4m_{eff}A_{loop}}{\hslash}\overset{.}{\phi}}},} & (24) \end{matrix}$ where N(t) is the number of times a polariton traverses the loop in time t. For the pSQUID ring geometry of radius r,

$\begin{matrix} {{{N(t)} = \frac{\hslash\; k_{0}t}{2\pi\;{rm}_{eff}}},} & (25) \end{matrix}$ where the momentum mode k₀=mv/

. Therefore, in a pSQUID the Sagnac phase is ϕ_(Sagnac)(t)=2k ₀ r{dot over (ϕ)}t,  (26) where r˜100 [μm] and k₀˜10 [μm⁻¹], such that k₀r˜10³. When the laboratory frame of reference rotates in time, the density profile, i.e., quantum interference of the condensate reads |ψ(x,t)|²+|ψ(x,t)|² cos[2k ₀ r(ϕ+{dot over (ϕ)}t)],  (27) allowing for quantum gyroscopic measurements to be taken from imaging the condensate rotating as a function of time. The SNR is given by the ratio of the Sagnac phase shift to the shot noise as a function of the number of photons luminescing from the condensate per second

$\begin{matrix} {{SNR} = {\frac{\phi_{Sagnac}(t)}{\sqrt{1\text{/}N}}.}} & (28) \end{matrix}$

The sensitivity {dot over (ϕ)}_(min) is found by setting SNR=1. For the vortex-antivortex superposition in the pSQUID geometry

$\begin{matrix} {{\overset{.}{\phi}}_{\min} = {\frac{1}{2k_{0}{rt}\sqrt{N}}.}} & (29) \end{matrix}$ Polariton densities are ˜10⁹ [cm⁻²], and occupy an area of ˜100 [μm²] such that N˜10³. As such, for the vortex-antivortex superposition in the pSQUID, {dot over (ϕ)}_(min)˜5×10⁻¹¹[rad·s ⁻¹·Hz^(−1/2)].  (30) From Eq. (13) the phase difference across the Josephson junctions is |ϕ₁(t)−ϕ₂(t)|=2πn−2k ₀ r{dot over (ϕ)}t.  (31) Moreover, using Eq. (14), the current amplitude in the pSQUID is I(t)=(I ₁ +I ₂)

  (32) where η=(I₁−I₂)/(I₁+I₂) describes the asymmetry across the Josephson junctions. Alternatively, Eq. (14) can be written using Eq. (23) as I(t)=(I ₁ +I ₂)

  (33) or Eq. (14) can be written in terms of the Gouy phase Eq. (17),

$\begin{matrix} {{I(t)} = {\left( {I_{1} + I_{2}} \right){\sqrt{{\cos^{2}\left\lceil {\frac{{\ell_{2}} - {\ell_{1}}}{2\left( {\ell_{2} - \ell_{1}} \right)}\frac{\overset{.}{z}t}{z_{R}}} \right\rceil} + {\eta^{2}\mspace{14mu}\sin^{2}\left\lceil {\frac{{\ell_{2}} - {\ell_{1}}}{2\left( {\ell_{2} - \ell_{1}} \right)}\frac{\overset{.}{z}t}{z_{R}}} \right\rceil}}.}}} & (34) \end{matrix}$

Allan variance is a method that is used to measure frequency stability in oscillators, or to determine the intrinsic noise in a system as a function of the averaging time. The Allan variance provides information related to five basic noise terms relevant to inertial sensor data. Namely, these noise terms are quantization noise, Angle Random Walk (ARW), Velocity Random Walk (VRW), Bias Instability (BI), rate random walk, and rate ramp. The Allan variance analysis of a time domain signal Ω(t) consists of computing its root Allan variance or Allan deviation as a function of different averaging times τ. Here, the time domain signal will be given by Eq. (22) for the accelerometer, and Eq. (30) for the gyroscope. The Allan variance can be computed in terms of the output angles θ, corresponding to each gyro, or accelerometer rate sample. As such,

$\begin{matrix} {\mspace{79mu}{{\theta(t)} = {\oint^{t}{{\Omega\left( t^{\prime} \right)}{{dt}^{\prime}.}}}}} & (35) \end{matrix}$

These angle measurements are made at times t=kτ₀ where k varies from 1 to N, i.e., (τ₀, 2τ₀, 3τ₀, . . . , Nτ₀). For a discrete set of samples, a cumulative sum can also be used to give N values of θ. Here the cumulative sum of the gyro or accelerometer output samples at each kτ₀ is taken and each sum obtained is then multiplied by sample period τ₀ to give N values of θ. Once N values of θ have been computed, the Allan variance is calculated using equation

$\begin{matrix} {{\sigma^{2}(\tau)} = {\frac{1}{2\tau^{2}}\left\langle \left( {\theta_{k + {2m}} - {2\theta_{k + m}} + \theta_{k}} \right)^{2} \right\rangle}} & (36) \end{matrix}$ where σ²(τ) represents Allan variance as a function of τ and

is the ensemble average. Upon expanding the ensemble average in Eq. (36),

$\begin{matrix} {{\theta^{2}(\tau)} = {\frac{1}{2{\tau^{2}\left( {N - {2m}} \right)}}{\sum\limits_{k = 1}^{N - {2m}}\;\left( {\theta_{k + {2m}} - {2\theta_{k + m}} + \theta_{k}} \right)^{2}}}} & (37) \end{matrix}$ where N is the total number of samples, m is the averaging factor, τ=mτ₀ is the averaging time and K is a set of discrete values varying from 1 to N. Eq. (37) computes the final rate Allan variance (e.g., via overlapping method) value using output angle θ for one particular value of τ. Once the output angle values of θ are known, the sample period and the value of m, Eq. (37) can be used to compute Allan variance. Finally, taking the square root of the result obtained from (37) to obtain the value of the root Allan variance or the Allan deviation for a particular value of τ. This result characterizes the noise in the device, AD(τ)=√{square root over (AVAR(τ))}.  (38)

White noise (i.e., random walk) appears on an Allan variance plot where the slope is −0.5. The random walk measurement (ARW for a rate-gyroscope, VRW for an accelerometer) is obtained at τ=1. From FIG. 7, it can be seen that for the polariton accelerometer VRW˜10⁻⁵[m·s ⁻¹·Hz^(−1/2)].  (39)

Allan variance curve (dots) in FIG. 7 correspond to the polariton accelerometer in the presence of AWGN with a SNR of 1 (i.e., 0 dB). The BI and VRW values are distinctly marked on the figure.

FIG. 8 shows the Allan variance curve (dots) corresponding to the polariton gyroscope in the presence of AWGN with a SNR of 1 (i.e., 0 dB). The BI and ARW values are distinctly marked on the figure.

Moreover, from FIG. 8, it can be seen that for the polariton gyroscope ARW˜10⁻⁵[°·s ⁻¹·Hz^(−1/2)].  (40)

BI appears on an Allan Variance plot at the minimum. For a pIMU, the BI measures how the bias of the polariton accelerometer or polariton gyroscope changes over a specified period of time at constant temperature. From FIG. 7, it can be seen that the BI for the polariton accelerometer is BI _(accel)˜10⁻³[m·hr⁻¹].  (41)

Moreover, from FIG. 8, it can be seen that the BI for the polariton gyroscope is BI _(gyro)˜10⁻³[°·hr⁻¹].  (42)

Results were obtained using the MATLAB code found in section VI, where the number of samples N=10⁴, and the sample period τ₀=1 sec.

TABLE I Ω_(min) BI ARW Sensitivity Technology [rad · s⁻¹ · Hz^(−1/2)] [° · hr⁻¹] [° · s⁻¹ · Hz^(−1/2)] Enhancement Polariton BECs   5 × 10⁻¹¹ 0.004385 1 × 10⁻⁵ 720,000 Ultracold Atoms   6 × 10⁻¹⁰ 0.04564 2 × 10⁻⁴ 72,000 Ring Laser Gyro (1 m²) 1.3 × 10⁻⁹ 0.1067 4 × 10⁻⁴ 6,000 Superfiuid Helium   2 × 10⁻⁷ 16.82 0.06535 4,000 Fiber Optic Gyro   3 × 10⁻⁷ 25.56 0.1004 26 (Commercial) Nanofabricated Gratings 3.6 × 10⁻⁶ 281.2 1.192 12 Moiré Classical Fringes 3.6 × 10⁻⁵ 2603 11.9 1

In classical electronic computing, information is usually encoded and accessed by adjusting and monitoring voltages on transistor switches inside integrated circuits. Each electronic transistor is set to a voltage state of either 0 or 1. In this way, electrical voltage is used to encode bits of information. As discussed above, polaritons are a neutral species, thereby disqualifying them for use as voltage state representations. As such, a charged supercurrent is needed to replace the voltage state representations of the classical transistor, as the supercurrent can represent a voltage state of 0 and 1 simultaneously, via superposition. In this way, room-temperature quatron superconductors can be used to build quantum transistors that replace the classical electronic transistors.

Quatrons are hybrid or quasi-particle optical-matter wave systems, which enable the quantum properties of light to be exhibited by matter waves. For example, as discussed above, OAM states are properties of light, which represents a Hilbert space. The Hilbert space found in the phase of photons allows quatron condensates to encode quantized phase information. Like polaritons, quatrons can exhibit BEC states at room-temperatures due their small effective masses. This allows matter-waves of room temperature superconducting quatrons.

Quatrons can therefore act as qudits, or quantum units of information. When the value of “d” is equal to two, qudits act as qubits, and can be made to represent and replace the binary voltage states of the classical transistor. Thus, quatron based quantum transistors are comparable to the classical electronic transistors. While classical transistors permit the encoding of two different voltage states one at a time, quantum transistors instead encode two or more magnetic field states simultaneously. This encoding of states is known as superposition and can be achieved in quatrons by using the photon coupling of OAM states as shown in connection with FIG. 2.

The superposition of the d-OAM states as magnetic fields allow the quatron OAM states to act as room-temperature quantum transistors. The magnetic field states correspond to the superposition states of the quatron condensate, originating from the OAM photon coupling. The quantum information can then be stored in a magnetic memory element local to each device. After the quantum computation is completed, the qubits can be measured and collapsed into their classical 0 or 1 states. The readouts can be used to query the value held by each qubit, and return the answer as a bit string of 0s and 1s.

Quatrons are composed of a hole and three-electrons coupled to photons in triple quantum layers inside semiconductor microcavities, for example. This makes qSQUIDs an ideal solution for replacing the classical electronic transistors, as quatrons are optically accessible bosons, capable of condensing at room-temperatures. These quatron qudits or qubits can be constructed into larger logical arrangements, and these larger logical arrangements can be made into larger computing elements that permit interaction, logical and mathematical operations, and the exchange and storage of information. In this way room-temperature programmable quantum computing can be achieved.

The quatron computing elements described are different from the classical Von Neuman architecture. One reason for this difference is that each qudit or qubit is associated with its own memory instead of centralizing such memory elements in, for example, a processor cache. In this way, the quatron based computing elements are neuromorphic in that the qudits or qubits behave more as neurons and the couplers accordingly act as synapses that control the flow of information between those neurons. This neuromorphic behavior will be beneficial to applications such as machine learning, as patterns can be instantiated in the architecture itself resulting in feedback loops between “hardware” and “software” in the devices.

FIG. 2 gives an example of an experimental configuration for the room-temperature qSQUID 200. An incident Gaussian beam of photons is passed through the waveguide to produce OAM vortex superpositions of light. This light is then passed through the Bragg Mirror configuration to produce a quatron-polariton Bose-Einstein condensate. Josephson junctions can be placed in a ring geometry inside the distributed Bragg reflector configuration, such that the device generates magnetic flux quantum Φ₀, capable of being measured by the Hall effect sensor 206 at the output. If the 2DEG component is removed from the system, the supercurrent will be charge neutral, thereby making it a polariton superfluid quantum interference device, in which case the Hall effect sensor 206 will need to be replaced with a photon detector, e.g., a CCD.

From FIG. 2 it can be seen that a Gaussian beam 202 is passed through a variable amplitude splitter 204 and phase shifter 203 such that the OAM emitter 201 generates an OAM superposition state of light. The OAM superposition state is then passed through a triple quantum layer 205 such that a quatron BEC is made. Owing to the phase of the OAM superposition state, the quatron BEC in the triple quantum layer 205 will obtain a vortex-antivortex superposition. Moreover, a circular channel can be placed in the triple quantum layer 205 such that a superconducting loop is obtained. Josephson junctions can be placed in the superconducting loop in the triple quantum layer 205 such that magnetic flux quanta are generated. The magnetic flux quanta can then be detected at the output of the device via the Hall effect sensor 206. Furthermore, the input coil and pickup loop 208 can be used to sense magnetic flux at the output of the device. In the quatron system, or the quatron superconducting quantum interference device (qSQUID) 200 as seen in FIG. 2, letting D be the distance between the layers of the triple-layered microcavity, r₁=r_(h)−r_(e1), r₂=r_(h)−r_(e2), and r₃=r_(h)=r_(e3), where r_(h) is the radius vector of the hole, and r_(e1) is the radius vector of the electron in the 2DEG layer, and r₂ and r₃ are the radius vectors of the electrons in the BEC layer, μ₁=m_(e1)+m_(h)/(m_(e1)+m_(h)), μ₂=m_(e2)m_(h)/(m_(e2))+m_(h), m_(h), m_(e1), and m_(e2) are the effective masses of the hole and electron in the 2DEG layer and the electron in the BEC layer, then the Hamiltonian reads

$\begin{matrix} {\hat{H} = \left\lceil \begin{matrix} {\hat{H}}_{coul} & {V\;{\delta\left( r_{1} \right)}} \\ {V\;{\delta\left( r_{1} \right)}} & {\hslash\omega}_{c} \end{matrix} \right\rceil} & (43) \end{matrix}$ where

ω_(c) is the energy of the bare cavity mode, and V describes the excitors-photon coupling. Furthermore,

$\begin{matrix} {{{\hat{H}}_{coul} = {{\sum\limits_{{i = 1},2,3}{\frac{\hslash^{2}}{2\mu_{i}}\frac{1}{r_{i}}\frac{\partial}{\partial r_{i}}\left( {r_{i}\frac{\partial}{\partial r_{i}}} \right)}} + {\frac{1}{r_{i}^{2}}\frac{\partial^{2}}{\partial\phi_{i}^{2}}} - \frac{e^{2}}{ɛ\; r_{1}} - \frac{e^{2}}{ɛ\sqrt{D^{2} + r_{2}^{2}}} - \frac{e^{2}}{ɛ\sqrt{D^{2} + r_{3}^{2}}} + \frac{e^{2}}{ɛ\sqrt{D^{2} + r_{1}^{2} + r_{2}^{2} - {2r_{1}r_{2}\mspace{14mu}{\cos\left( {r_{1},r_{2}} \right)}}}} + \frac{e^{2}}{ɛ\sqrt{D^{2} + r_{1}^{2} + r_{3}^{2} - {2r_{1}r_{3}\mspace{14mu}{\cos\left( {r_{1},r_{3}} \right)}}}} + \frac{e^{2}}{ɛ\sqrt{\left( {2D} \right)^{2} + r_{3}^{2} + r_{2}^{2} - {2r_{3}r_{2}\mspace{14mu}{\cos\left( {r_{3},r_{2}} \right)}}}}}},} & (44) \end{matrix}$ where e is the elementary charge. The vector for the superposition of OAM states in the quatron system is of the form

$\begin{matrix} {{\left\lceil \begin{matrix} {{\alpha\psi}\left( {r_{i},\phi_{i},t} \right)} \\ \beta \end{matrix} \right\rceil = \left\lceil \begin{matrix} {\alpha\frac{\left. \left. {{\left( {❘\ell} \right\rangle +}❘{- \ell}} \right\rangle \right)}{\sqrt{2}}} \\ \beta \end{matrix} \right\rceil},} & (45) \end{matrix}$ where α=−β/√{square root over (S)}=√{square root over (2)}/2, S is the surface area, and |

is given by Eq. (1). A charged supercurrent device can be obtained by implementing the Josephson junctions with quatrons. For the case when the quatron system, i.e., the qSQUID as shown in FIG. 2 contains two Josephson junctions, the superconducting phase evolution is given by

$\begin{matrix} {{{U(t)} = {{- \frac{\hslash}{2e}}{{{\phi_{1}(t)} - {\phi_{2}(t)}}}}},} & (46) \end{matrix}$ owing to the quatron charge −2e. Rewriting Eq. (46) using Eq. (31) for the ring geometry, such that

$\begin{matrix} {{U_{k_{0}}(t)} = {{- \frac{\hslash}{e}}\left\lceil {{\pi\; n} - {k_{0}r\;\overset{.}{\phi}\; t}} \right\rceil}} & (47) \end{matrix}$ where n is integer valued, and represents the periodicity of the superconducting loop, as in Eq. (13). Alternatively, in terms of the (superconducting) magnetic flux quantum Φ₀, rewriting Eq. (47) as

$\begin{matrix} {{U_{k_{0}}(t)} = {{- \Phi_{0}}{\left\lceil {n - \frac{k_{0}r\;\overset{.}{\phi}\; t}{\pi}} \right\rceil.}}} & (48) \end{matrix}$

This allows the measurement of the superconducting phase as a function of the Sagnac phase in the ring geometry.

A qubit (i.e., a two-level quantum system), evolves in a Hilbert space spanned by two basis vectors usually denoted |0

and |1

in analogy with classical information. For a pare state, it can be represented by: |ψ

=a|0

+b|1

  (49)

Where a and b are complex: and |a|²+|b|²=1. In terms of polariton BECs created via OAM states of light carrying quantized OAM

, defining the corresponding logical basis vectors as |0

=|+

=|R

and |1

=|−

=|L

. Here, R and L refer to the right and left handedness of the helical wavefront, respectively. In analogy to the case of a polarization basis, the specific following modes as vertical, diagonal, horizontal and anti-diagonal are denoted as:

$\begin{matrix} {{\left. {{\left. {❘V} \right\rangle = \frac{\left. {{\left. {❘R} \right\rangle -}❘L} \right\rangle}{\sqrt{2}}},{❘D}} \right\rangle = \frac{\left. {{\left. {❘R} \right\rangle + i}❘L} \right\rangle}{\sqrt{2}}},} & (50) \\ {{\left. {{\left. {❘A} \right\rangle = \frac{\left. {{\left. {❘R} \right\rangle - i}❘L} \right\rangle}{\sqrt{2}}},{❘H}} \right\rangle = \frac{\left. {{\left. {❘R} \right\rangle +}❘L} \right\rangle}{\sqrt{2}}},} & (51) \end{matrix}$ The bases |R

, |L

, |H

, |V

, and |D

, |A

will constitute the 3 mutually unbiased bases necessary to form a Bloch sphere. If the normalization is written implicitly and the global phase is ignored, Eq. (49) can now be rewritten as

$\begin{matrix} \left. {{\left. {\left. {❘\Psi} \right\rangle = {{\cos\left( \frac{\theta}{2} \right)}❘R}} \right\rangle + {{\sin\left( \frac{\theta}{2} \right)}{\exp\left( {i\;\phi} \right)}}}❘L} \right\rangle & (52) \end{matrix}$ For example, in a Franson interferometer configuration, an entangled state can be obtained, such that

$\begin{matrix} {\left. {❘\Phi} \right\rangle = {\frac{\left. \left. {\left. {{\left. {\left. {❘L} \right\rangle_{A}❘L} \right\rangle_{B} + {\exp\left( {\phi_{A} + \phi_{B}} \right)}}❘R} \right\rangle_{A}❘R} \right\rangle_{B} \right)}{\sqrt{2}}.}} & (53) \end{matrix}$ Once the charged quatron superconductor inherits an OAM

, a quantized flux will be generated. As such, the superconducting phase can be written in terms of the OAM,

$\begin{matrix} {{U_{\ell}(t)} = {{- \Phi_{0}}\left\lceil {n - {\ell\frac{\overset{.}{\phi}t}{\pi}}} \right\rceil}} & (54) \end{matrix}$ to induce a flux quantum, where

is the OAM winding number, {dot over (ϕ)} is the angular velocity of the condensate, and n is integer valued, representing the periodicity. In this context, it is useful to discuss the various types of qubits (i.e., quantum bits). There are two primary types of qubits that exploit the superconductivity of Cooper pairs (i.e., electron-electron pairs), namely, charge qubits and flux (i.e., phase) qubits. The Josephson junction is the basic building block of a superconducting qubit, and consequently the basic entity of a quantum computer or a quantum simulator. The Josephson junction consists of two superconducting regions separated by a “weak link.” Following this reasoning, a quatron superconducting flux qubit, where a pair of conductivity electrons is linked with an exciton-polariton by Coulomb interaction, and used analogously to the Cooper pairs of the superconducting flux qubit. In this case, a quatron flux qubit consists of a ring of one or more Josephson junctions. For the case of a ring of superconducting quatron BEC, two Josephson junctions are implemented. For example, using orbital angular momentum of light, one can encode |0

and |1

data into these counter-rotating quatron superconducting states, effectively constituting a novel quatron flux qubit using Eq. (16). Herein, the advantages of room-temperature quatron BEC and Josephson junction based implementations are combined. The basic idea is to use a ring of quatron BEC to realize a room-temperature quatron BEC analogue of the superconducting flux qubit and similar technologies. Supercurrents flowing in a ring provide a physical implementation of a qubit, as they generate a flux characteristic to the state of the qubit being measured. In this way, quatron flux qubits with supercurrents flowing in a ring-shaped geometry rotated by Josephson junction(s) induce a flux. In quatron condensates, the supercurrent is caused by a gradient in the macroscopic phase

$\begin{matrix} {{J_{s} = {{- 2}e\;\rho_{s}\left\lceil {{\frac{\hslash}{m_{eff}}{\nabla{\phi_{ext}\left( {x,t} \right)}}} + {\frac{2e}{m_{eff}}{A\left( {x,t} \right)}}} \right\rceil}},} & (55) \end{matrix}$ where A(x,t) is the magnetic vector potential, and m_(eff) is the quatron effective mass. Furthermore, in a circular loop, Eq. (55) can only change as a multiple of 2π, meaning that the supercurrent is always quantized in this particular geometry, as seen in FIG. 6. Flux quantization is the direct result of the macroscopic coherence that is inherent to BECs. Condensates of quatron-polaritons exist at room temperatures. A quatron is a pair of conductivity electrons coupled with an exciton-polariton via Coulomb interaction, and thus a quatron carries charge −2e. Because quatrons are optically accessible bosons, they can be coupled with the photon OAM to create states as seen in Eq. (49).

INS generally involve using rotation (gyroscope), and acceleration (accelerometer) sensors to determine a location relative to a given starting point. Some of the applications for INS include, but are not limited to attitude and heading reference systems (AHRS), Global Position Systems (GPS), submarines, aerospace vehicles, petroleum exploration, drilling, and automobiles. An INS typically requires three accelerometers and three gyroscopes one for each physical axis. The accelerometers can be used to determine tilt based on the projection of gravitational acceleration onto each axis.

The GPS satellite system is not sensitive to angular orientation or rotation rate and so gyroscopes are still required for many applications, even where GPS can be utilized. Differential GPS utilizes additional information supplied by local transmitters and therefore can be used to enhance accuracy. Moreover, and particularly for military applications, the signals from GPS satellites are thought to be vulnerable to jamming and the satellites are themselves vulnerable to attack. Additionally, satellite signals cannot be received under water, dense urban environments, or even where heavy-tree cover is encountered. Accordingly, there remains a need for sensitive INS systems that do not rely upon the GPS satellite system, if even only as a supplement. In such conditions, even small biases in the readings of accelerometers or gyroscopes can result, over time, in increasing errors in navigation.

FIG. 9 depicts a strapdown IMU (sometimes also referred to as an INS) navigation algorithm integrating the polariton accelerometers and polariton gyroscopes as discussed above into a pIMU 800.

Inertial navigation is a technique whereby measurements provided by accelerometers and gyroscopes are used to track the position and orientation of an object relative to a known starting point, orientation and velocity. These know starting points, orientations, and velocities typically represent the results of the calculations made at an initial time. Inertial navigation units typically include three at least partially orthogonal gyroscopes 801 and three at least partially orthogonal accelerometers 802 that measure angular velocity and linear acceleration respectively. By processing the measurements obtained from these devices it is possible to track the position and orientation of an object in three dimensional space.

To track the orientation of an object, the measurements obtained from the three axes of at least partially orthogonal gyroscopes 801 are integrated 803 using well-known mathematical techniques. To track the position of the object, three at least partially orthogonal accelerometers 802 are used. The measurements obtained from the three at least partially orthogonal accelerometers 802 are resolved into global coordinates or global axes 804 using the known orientation as determined by the integration of the gyroscope 801 measurements, again using well-known mathematical techniques.

Next, the IMU 800 corrects for the force of gravity on the object in block 805. The output of this set of calculations is the global acceleration of the object in three dimensional space. This global acceleration is integrated in block 806 using the initial velocity (calculated at the previous time interval) to obtain the current velocity. Next, the current velocity is integrated with the initial position (again, calculated at the previous time interval) to obtain the current position. The values of orientation, velocity and position are typically stored to be used for the calculation of orientation, velocity and position for the next time interval. Alternatively, the higher order derivatives of linear velocity, angular velocity and position can be taken.

As can be appreciated by those of ordinary skill in the art, the three at least partially orthogonal gyroscopes 801 can be constructed using the polariton gyroscopes discussed above. Likewise, the three at least partially orthogonal accelerometers 802 can be constructed using the various forms of the polariton accelerometers discussed above.

The resulting polariton-based gyroscopes and accelerometers of the present disclosure achieve high degrees of accuracy and sensitivity. Allan variance is a method that is used to measure the frequency stability in oscillators, or to determine the intrinsic noise in a system as a function of the averaging of time. The Allan variance provides information related to five basic noise terms relevant to inertial sensor data. Namely, these noise terms are quantization noise, angle random walk, bias instability, rate random walk, and rate ramp. In addition to navigational applications, the accelerometers and/or gyroscopes of the present disclosure can be used for other important applications that require exquisite sensitivity, such as determining small changes in the rotation rate of the Earth and the associated length of day fluctuations. Such fluctuations are caused primarily by tidal forces that alter the mass distribution of the oceans that cover most of the surface of the Earth. As a result, the rotation rate of the Earth (and thus the length of the day) alters correspondingly to conserve angular momentum. Likewise, the devices of the present disclosure can also be used in seismic monitoring. Little information exists regarding the rotational spectrum of earthquakes and other seismic events. Sensitive measuring devices such as described in the present disclosure can provide greater depth of information over currently used techniques, such as Very-long-baseline interferometry (VLBI) because they are much more portable than VLBIs and therefore can take local measurements.

Another application is in the field of torsion pendulums which are used, for example, to take measurements of the gravitational constant G. The devices described in the present disclosure are sensitive enough to allow the detection of thermal noise (that results in the geometry of the pendulum itself) and reduce the resulting systematic error.

Finally, the accelerometers and gyroscopes of the present disclosure have application in the field of general relativity testing the extent and limits of the Lense Thirring effect and the impact of geodetic precession. Lense and Thirring predicted in 1918 that a satellite orbiting a rotating body would experience a precession of its orbit due to the rotating mass. The Lense-Thirring effect can be obtained in the limit of weak gravitational fields by considering the relativistic field equations in the linear approximation, and leads to off-diagonal terms in the metric. The effect on nearby space-time is similar to the behavior of a viscous fluid surrounding a rotating sphere, and is sometimes called the dragging of inertial frames.

It is believed the examples of the present disclosure and their advantages will be understood from the foregoing description, and it will be apparent that various changes may be made thereto without departing from the spirit and scope of the disclosure or sacrificing all of its advantages, the examples hereinbefore described merely being preferred or exemplary.

Those of ordinary skill in the art will appreciate any of the aforementioned techniques may be used without departing from the present disclosure.

VI. Appendix 1 % Allen Deviation 2 function [T, sigma] = allan (omega, fs, pts) 3 [N,M] = size (omega); % figure out how big the output data set is 4 n = 2. ^(∧)(0:floor(log2(N/2)))′; % determine largest bin size 5 maxN = n(end); 6 endLogInc = log10 (maxN); 7 m = unique (ceil (logspace (0, endLogInc, pts)))′; % create log spaced vector average factor 8 t0 = 1/fs; % t0 = sample interval 9 T = m★t0; % T = length of time for each cluster 10 theta = cumsum(omega)/fs; % integration of samples over time to obtain output angle 11 sigma2 = zeros (length (T), M); % array of dimensions (cluster periods) X (# variables) 12 for i=1:length(m) % loop over the various cluster sizes 13  for k=1:N−2★m(i) % implements the summation in the AV equation 14  sigma2(i, :) = sigma2(i, :) + (theta(k+2★m(i), :) − 2★theta (k+m (i), :) + theta(k, :)). ^(∧)2; 15  end 16 end 17 sigma2 = sigma2./repmat ((2★T. ^(∧)2. ★ (N−2★m)) , 1, M); 18 sigma = sqrt(sigma2); 1 % Accelerometer 2 t=1:10^(∧)4; 3 % additive white Gauesian noise 4 omega=transpose(awgn (5★10 ^(∧) (−9) ★t, 0, ′measured′)); 5 [T, sigma] = allan (omega, 1, 10^(∧)4); 6 loglog (T, sigma); 7 find(sigma==min(sigma)) % find BI point 1 % Gyroscope 2 t=1:10^(∧)4; 3 % additive white Gauesian noise 4 omega=transpose (awgn ((180/pi) ★5★10^(∧) (−11). ★t,0, ′measured′)); 5 [T, sigma] = allan (omega, 1, 10^(∧)4); 6 loglog (T,sigma); 7 find (sigma==min (sigma)) % find BI point 

What is claimed is:
 1. A quasi-particle superfluid Inertial Measurement Unit (IMU) comprising: one or more accelerometers, using quasi-particle Bose-Einstein condensation as a lasing source, arranged at least partially orthogonally; one or more gyroscopes using quasi-particle Bose-Einstein condensation as a lasing source; arranged at least partially orthogonally; one or more processors operable to execute instructions stored in a memory, the one or more processor operable to: determine an orientation of the quasi-particle superfluid IMU based on the outputs of the one or more gyroscopes, wherein the orientation is determined by integrating, using an IMU algorithm, the outputs of the one or more gyroscopes; determine a velocity of the quasi-particle superfluid IMU based on the outputs of the one or more accelerometers; determine a position of the quasi-particle superfluid IMU from the orientation and velocity based on previously known velocity and position or higher order derivatives of linear velocity, angular velocity, and/or position, wherein the outputs of the one or more accelerometers are resolved into coordinates or axes using the orientation of quasi-particle superfluid IMU.
 2. The quasi-particle superfluid IMU of claim 1, wherein the quasi-particle superfluid is comprised of exciton-polariton Bose-Einstein condensation.
 3. The quasi-particle superfluid IMU of claim 1, wherein at least one of the one or more accelerometers comprises a Mach-Zehnder interferometer and/or Franson interferometer configuration.
 4. The quasi-particle superfluid IMU of claim 1, wherein at least one of the one or more gyroscopes comprises a Mach Zehnder interferometer and/or Franson interferometer configuration.
 5. The quasi-particle superfluid IMU of claim 1, wherein at least one of the one or more accelerometers or at least one of the one or more gyroscopes comprise an interferometer further comprising a device selected from the group consisting of spiral phase plates, spatial light modulators, cylindrical lens mode converters, and/or computer generated holograms.
 6. The quasi-particle superfluid IMU of claim 1, wherein at least one of the one or more accelerometers or at least one of the one or more gyroscopes comprise a Dove prism.
 7. The quasi-particle superfluid IMU of claim 1, wherein at least one of the one or more accelerometers or at least one of the one or more gyroscopes comprise a pumped microcavity quantum well device.
 8. The quasi-particle superfluid IMU of claim 1, wherein a pumped microcavity quantum well device is selected from the group consisting of photonic crystals, optical microcavities, planar dielectric Fabray-Perot optical cavities, polymers, triple layer structure, and semiconductor microcavities.
 9. The quasi-particle superfluid IMU of claim 1, wherein a Bose-Einstein condensate is measured by a detector to determine a property of the Bose-Einstein condensate selected from the group consisting of the phase, the flux quanta, the linear velocity or the angular velocity.
 10. The quasi-particle superfluid IMU of claim 1, wherein a Bose-Einstein condensate is quantum entangled with superposition states of light.
 11. The quasi-particle superfluid IMU of claim 10, wherein the superposition states of light are squeezed superposition states of light.
 12. The quasi-particle superfluid IMU of claim 10, wherein the superposition states of light are vortex-antivortex superposition states of light.
 13. The quasi-particle superfluid IMU of claim 1, wherein the quasi-particle superfluid IMU is operable at room temperature.
 14. A method of a quasi-particle superfluid Inertial Measurement Unit (IMU), the method comprising: determining an orientation of the quasi-particle superfluid IMU based on an output of one or more gyroscopes, wherein the quasi-particle Bose-Einstein condensation is used as a lasing source, arranged at least partially orthogonally, wherein the orientation is determined by integrating, using an IMU algorithm, the outputs of the one or more gyroscopes; determining a velocity of the quasi-particle superfluid IMU based on an output of one or more accelerometers, wherein the quasi-particle Bose-Einstein condensation is used as a lasing source, arranged at least partially orthogonally; determining a position of the quasi-particle superfluid IMU from the orientation and velocity based on a previously known velocity and position or higher order derivatives of linear velocity, angular velocity, and/or position, wherein the outputs of the one or more accelerometers are resolved into coordinates or axes using the orientation of quasi-particle superfluid IMU.
 15. The method of claim 14, wherein the quasi-particle superfluid is comprised of exciton-polariton Bose-Einstein condensation.
 16. The method of claim 14, wherein the quasi-particle superfluid Inertial Measurement Unit (IMU) is operable at room temperatures.
 17. The method of claim 14, wherein a quasi-particle Bose-Einstein condensate is quantum entangled with superposition states of light.
 18. The method of claim 15, wherein the superposition states of light are squeezed.
 19. A quasi-particle superfluid Inertial Measurement Unit (IMU) comprising: one or more quasi-particle Bose-Einstein condensation based accelerometer arranged at least partially orthogonally for determining the velocity of the superfluid IMU; one or more quasi-particle Bose-Einstein condensation based gyroscopes arranged at least partially orthogonally for providing the orientation of the superfluid IMU, wherein the orientation is determined by integrating, using an IMU algorithm, the outputs of the one or more gyroscopes; and one or more processor for determining from the orientation and velocity based on previously known velocity and position or higher order derivatives of linear velocity, angular velocity, and/or position, wherein the outputs of the one or more accelerometers are resolved into coordinates or axes using the orientation of quasi-particle superfluid IMU. 